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ABSTRACT 



• A computer program for simulating the performance of 
the 100-B Surface Effect Ship is used to study the longi- 
tudinal motions of the ship under various wave conditions. 
An investigation into the effect that waves have on bow 
seal, stern seal, and plenum presures is conducted. The 
relationship betvjeen the different pressures and their 
associated requirements of input air from the fans is 
studied. It is concluded that the computer program is 
limited to a certain range of speeds. The concept of the 
ship capturing air to reduce drag and increase its speed 
is questionable due to the rapid replenishment of air 
required to keep the ship riding on its bubble of air. 
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PSP 


= 


pounds per square foot 


PSI=1^ 


= 


yaw angle 


Q 


= 


rotational velocity about y-axis 


*^ln 


= 


air flow rate in 


^out 


r= 


air flov; rate out 


R 


= 


rotational velocity about z-axis 


RHS 


= 


right hand side 


SAM 


= 


shear and moment 


SIDV/L 


= 


sidewall 


STNSL 


= 


stern seal 




= 


standard atmospheric reference 


T 


= 


arbitrary steady-state time 


THETA=e 


= 


roll angle 


U 


= 


longitudinal velocity 


V 




lateral velocity 


W 




vertical velocity 


X 


= 


horizontal distance in direction of motion 


Y 


= 


horizontal distance perpendicular to direction 
of motion 


z 


zz 


vertical distance 
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I . INTRODUCTION 



Conventional monohull displacement ships are limited to 
a maximum speed of approximately 50 knots. This is because 
the drag on conventional displacement ships Increase exponen- 
tially with speed, resulting in a limit as to maximum 
economical propulsion plant that can be installed. 

Many new designs of ships have been proposed over the 
past few years. These Include the "slender body" ship, 
catamaran, hydrofoils, semlsubmerged ships, air cushion 
vehicles (ACV), and surface effects ships (SES). All of 
these types of ships or vehicles have certain advantages 
and disadvantages. 

One of the most promising designs in overcoming the 50 
knot speed limit is the surface effect ship. The surface 
effect ships are usually classified as to either the Captured 
Air Bubble (CAB) type or the Hovercraft (GEM-Ground Effects 
Machine) type. 

The Hovercraft design is size limited due to the 
Inefficient use of forced air to keep the craft above the 
water surface. No attempt is made to "capture" the air 
under the craft although the design of flexible skirts 
around the craft have been improved in recent years to 
retard the leakage rate of the forced air. 

The Captured Air Bubble design uses rigid sidewalls and 
flexible seals (bov; and stern) to "capture" the forced air 
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in the plenum, thereby reducing the size of the machinery 
required to generate the input air flow. The drag created 
by the CAB craft is reduced considerably when the craft is 
operating in the "on bubble" mode. This in turn allows for 
a greater speed than conventional monohull ships for the 
same amount of thrust. The CAB is not size limited as is 
the case of the Hovercraft. 

To date three CAB type vehicles are operational. 

a. The XR-3 craft was built by David Taylor Model 
Basin (now designated Naval Ships Research and Development 
Center [NSRDC]). It has a displacement of approximately 
2.5 tons. The XR-3 is presently located at the Naval 
Postgraduate School (NFS), Monterey, California. 

b. The SES 100-A craft v;as built by Aerojet-General 
Corporation. It has a displacement of approximately 100 
tons and is located at Tacoma, Washington. 

c. The SES 100-B craft was built by Bell Aerospace 
Company with headquarters at Michoud Station, Louisiana. 

It has a displacement of approximately 100 tons and is 
presently located at Lake Pont chart rain, Louisiana. To 
date this craft has been run at speeds in excess of 70 
knots [Ref. 1]. 

To better understand the motions and loads of the SES 
100-B CAB, the Surface Effect Ships Project Office in 
Washington, D.C. issued a contract to Oceanics, Inc. of 
New York to write a computer program and gather certain 
basic data on the SES 100-B. The computer program was 
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completed In 1971 and a copy of the program v;as Installed 
in the W.R. Church Computer Center at Monterey in October, 
1972. 

This investigation will concern Itself only v;ith pre- 
senting results obtained in exercising the computer program 
developed for the SES 100-B CAB. 

In a letter to Doctor G.J. Thaler, NFS, Monterey, from 
Mr. S. Davis, SESPO, [Ref. 2], a list of future studies of 
the SES 100-B craft, using the computer program developed 
by Oceanics, was suggested. This author selected to study 
the variation of pitch, heave, and pressures in regular 
waves (ahead and astern) by varying ship speed, wave length, 
and wave height as suggested in paragraph 2a, [Ref. 2]. 
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II. NATURE OP THE PROBLEM 



A. LEADING PARTICULARS 

The following is a list of the leading particulars of 
the SES 100-B craft. 

1. Dimensions 



Length, Overall 
Beam, Overall 
Height, Overall 
Hull Length 
Cushion Area 
Bubble Pressure 

Distance Centerline to Sidewall 
Plenum Length at Water Surface 
Plenum V/idth at Water Surface 
Plenum Average Height 
Froude Number 
Keel to Deck 

2 . Weights 
Empty Weight 
Gross Weight 

3. Center of Gravity 

Longitudinal Feet Forward of Transom 
Transverse Feet to Starboard 
Vertical Feet above Keel 



76.3 ft. 

35.0 ft. 

27.0 ft. 

72.0 ft. 
1907 ft.' 

92.8 PSP 
15.5^' ft . 

65.31 ft. 
31.16 ft. 
6.17 ft. 
0.6 

13.5 ft. 



103,071 lbs 
209,999 lbs (105 tons) 



33.09 ft. 

0.0 f t . 

7.09 ft. 
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^ . Propulsion 



Three (3) Pratt and Whitney FT 12A-6 Marine Gas 
Turbines. il500 HP (each) Rated at 59®P. 

5 . Propellers 

Two (2) Supercavltating CRP Marine Screws, 



Diameter 


3.50 


ft. 


Feet Forward of Transom 


-1.0 


ft. 


Feet from Centerline to Propeller Center 


• 15.75 


ft. 


Feet Above Keel 


.09 


ft . 


Fans 

Number of Bow Seal Fans 


1 




Speed of Bow Seal Fans 


1700 


RPM 


Number of Cushion Fans 


5 




Speed of Cushion Fans 


1700 


RPM 


Number of Stern Seal Fans 


2 




Speed- of Stern Seal Fans 


1870 


RPM 


Maximum Total Fan Power 


1187.5 


HP 


Rudders 

Number of Rudders 


2 




Feet Forward of Transom 


6.125 


> ft 


Feet from Centerline to Rudder Centroid 


15.6 


ft 


Feet Above Keel 


-1.5 


ft 


Rudder Area (one) 


6.375 


> ft 


Appendages (fins) 
Number of Appendages 


2 




Appendage Span 


3.0 


ft 


Appendage Chord 


3.0 


ft 
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Distance of Appendage Centroid Forv;ard 



of Transom 



15.0 ft. 



Distance from Centerline to Centroid 



of Appendage 



1'I.875 ft. 



Distance of Appendage Centroid Above 



Keel 



-1.5 f t . 



Appendage Cant Angle (positive cant 
angle is inboard) 



30.0 deg. 



B. COORDINATE SYSTEM 

A body moving in a fluid can move in all six degrees of 
freedom of motion — translation along each of three 
orthogonal axes and rotation about each of the three axes. 

The SES 100-B craft has a plane of symmetry (as have most 
ships in existence), which is defined by the ship's center- 
line and a line perpendicular to the ship's deck. The axis 
system describing the SES 100-B craft is assumed to trans- 
late with the craft and only to allow rotation of the axes 
about the vertical axis. With the inertial frame axes being 
a right handed cartesian coordinate system, represented by 
X , y , and z , with z positive down, and an origin selected 
at a reference point on the undisturbed free surface, the 
z-axes are identical. The ship axis system has its origin 
at the center of gravity of the ship with the x-axis positive 
toward the bow, y-axis positive to starboard, and z-axls 
positive downvrard. The relationship between the inertial 
reference frame and craft axis system is based upon the 



values of pitch (6), roll (<!)), and yaw (H') angles. This 
relationship can then be expressed [Ref. 3] 






cos 


Y 


sin 


Y 


0 




X 


= 


-sin 


Y 


cos 


Y 


0 




Y 




0 




0 




1 




Z 



( 1 ) 



C. WAVE AMPLITUDE 

The computer program does not allow for an instantaneous 
build-up of wave amplitude. Instead an exponential rise 
in wave amplitude is used based on the formula 

- time 

wave amplitude = 1.0 - e 7 ^2) 



A one foot amplitude will require approximately ^.0 
seconds to build up to 95 percent of maximum amplitude and 
approximately 6.0 seconds to build up to 99 percent of 
maximum amplitude.. The reason for this exponential build-up 
of wave amplitude is that the integration routine used in 
the computer program cannot calculate accurately such an 
instantaneous Jump. 

D. COMPUTER SUBPROGRAMS 

The computer program developed for the SES 100-B craft 
uses various subprograms for calculating data. These 
subroutines are [Ref. 4]: 
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1 . Main Program 



The main program contains the logic for inter- 
connecting the various subroutines. It compares the running 
value of time with the finish time, calculates the next 
value of time for printing after each print time, and 
provides some output data. It contains the logic for 
integration of the trajectory (x,y,'i') by the trapezoidal 
rule . 

Output is time, U, V, W, P, Q, R, Z, Phi, Theta, 

X, Y, Psi, sideslip angle, rudder angle, and thrust. 

2 . Subroutine AEROD 

Subroutine AEROD calculates the aerodynamic 
forces and moments on the craft. 

Output is aerodynamic forces and moments. 

3 . Subroutine BOWSL 

Subroutine BOWSL calculates the forces and moments 
due to the bow seal and the leakage flow rate associated 
with any gaps under the various bow seal stations. 

Output is the gap at the various stations, wetted 
lengths at the various stations, and forces and moments due 
to the bow seal. 

^ . Subroutine COLFIL 

Subroutine COLFIL gives a tabular summary of 
selected variables in the vertical and/or lateral plane. 

A total of eight variables can be selected for each 
summary. 
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5 . Subroutine DMINV 



Subroutine DMINV Is a matrix Inversion package used 
for Inverting the mass matrix. The method used Is the 
standard Gauss-Jordan technique. The determinant of the 
matrix Is also calculated. 

There Is no output from this subroutine. 

6 . Subroutine FAN 

Subroutine PAN Is used to calculate the Inflov; to 
the plenum from the fans for a given pressure differential 
across the fans. 

Output Is flow rates and pressure differentials 
across the bow seal, plenum, and stern seals. 

7 . Function FGl 

The purpose of PGl Is to perform the task of 
evaluating a tabular function y = f(x). The procedure 
used Is linear Interpolation between data points. 

There Is no output data from function FGl. 

8 . Subroutine PORIT 

Subroutine PORIT Is used to obtain a Fourier analysis 
of a periodically tabulated function. 

There Is no output data from subroutine FORIT. 

9 . Subroutine INCON 

Subroutine INCON contains the logic for the reading 
of all Input data, the Initialization of variables, and 
for the Initiation of new cases. 
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Output is a complete listing of all initial condi- 
tions. The initial conditions will always be printed at 
the start of each new case. 

A sample data deck for establishing initial 
conditions is contained in Appendix B. 

10 . Subroutine INTGRL 

Subroutine INTGRL is used to integrate a system of 

first order ordinary differential equations. It uses a 

variable time step technique based on the Runge-Kutta-Merson 

algorithm. The program will stop the calculation if the 

—6 

time step becomes smaller than 10 sec. 

Outputs will occur when the calculated error is 
greater than the tolerance set on that particular integra- 
tor. Output format contains time, integrator step size 
(BELT), integrator number (J), error value, and tolerance 
value. 

11 . Subroutine PROP 

Subroutine PROP calculates the forces and moments 
on the craft due to the propulsion system. The port engine 
is always assumed ,to be the engine which fails if an engine 
out is designated. 

Output is the propulsion forces and moments. 

12 . Subroutine RHS 

Subroutine RHS contains the FORTRAN expressions for 
the right hand side of the system of first order differential 
equations. 
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Output consists of plenum pressure, fan power re- 
quired, fan flow rate, leakage flow rates for bov 7 , sidewall, 
and stern seals, plenum area and volume, total forces and 
moments, bov; and stern accelerations, and an array of 
values of the right hand side of the differential equations. 

13 . Subroutine RUDDER 

Subroutine RUDDER calculates the forces and moments 
on the rudder as well as the rudder motions. 

Output is the forces and moments on the rudder. 

1 ^ . Subroutine SAM 

Subroutine SAM performs all the processing and 
calculations involving determination of shears and moment. 

This subroutine was not required to be used in 
obtaining the data contained in this report. 

15 • Subroutine SIDEV/L 

Subroutine SIDEV/L calculates the forces and moments 
due to the sidev;alls as well as the leakage flow rates 
associated with any gaps under the sidewalls. The subroutine 
also contains the logic for the forces and moments due to 
the appendages (fins). 

Output is the gap and Immersion depth at designated 
stations along port and starboard sidewalls and also the 
forces and moments due to the sidewalls. 

16 . Subroutine STNSL 

Subroutine STNSL calculates the forces and moments 
acting on the craft due to the stern seal as well as the 
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leakage flov; rates arising from any gaps which open under 
the seal. 

Output is gap and stern seal v/etted lengths at 
various stations across the stern seal in addition to forces 
and moments arising from the stern seal. 

17 . Functions T1 and T2 

Functions T1 and T2 are two function subroutines 
used to calculate various trigonometric relations used 
by subroutine WAVES. 

There is no output data from Functions T1 and T2 . 

18 . Subroutine VfAVES 

Subroutine WAVES calculates the vrave forces and 
moments acting on the craft. It also generates the wave 
amplitudes at various stations around the seals and sidewalls, 
as well as bubble volume lost due to wave elevation. The 
sea state for irregular seas is computed by adding together 
a series of regular v;aves with an appropriate distribution 
of amplitude and frequency. 

Output is the' wave elevations relative to calm water 
at the various stations around the bow seal, stern seal, 
port and starboard sidewalls. Also the wave elevation at 
the center of gravity, plenum volume lost due to waves, 
and wave forces and moments are listed. 

E. EQUATIONS 

The ten (10) first order ordinary differential equations 
derived [Ref. 3] for use by a digital computer are; 
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m§: 











u 



mx 



(3) 



V = my mur 


(^) 


• 

v; = mz 


(5) 


p = I K - I 
^ XX xz 


(6) 


• 

q = I M 

yy 


(7) 


r = -I K + I 

xz zz 


(8) 


• 

w = z 


(9) 


• 

p = (j) 


(10) 


« 

q = e 


(11) 


• 

m, = p„(Q. - Q , ) 

b a in out 


(12) 



In addition three (3) auxiliary equations [Ref. 3] are 
used. These are 



t 

x=/udt (13) 

o 

t 

y = / V dt (1^) 

o 

t 

¥ = / r dt (15) 

o 



Since the variables x, y, and Y defined in Equations 
(13) - (15) are expected to have (in general) slow variation 
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and do often extend to large values, Simpson’s Rule 
Technique is used for integration in the Main Program 
section. 

Equations ( 3 ) - (12) use the Runge-Kutta-Merson method 
of integration. A tolerance level is selected for each of 
these equations. A truncation error, e, is calculated. If 
any element in e > tolerance level, then the trial time 
step size previously selected is halved and the computations 
are repeated for that time step. If all elements of 
e £ l/l6 tolerance, the step size is then doubled in the 
next integration step. Otherwise the computations proceed 
using the same step size for the next calculations. 

F. NUi^ERICAL INTEGRATION TECHNIQUE 

There are many integration schemes to choose from for 
a general system simulation v;here speed and accuracy is of 
paramount importance. The Runge-Kutta-Merson method was 
selected [Ref. 3]. ' 
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III. COMPUTER PROGRAM PROCEDURE 



A. PROGRAM FAMILIARITY 

The first priority in using the computer program devel- 
oped by Oceanics, Inc. was to become familiar with the 
various subroutines and to determine the limitations of 
the program. Very little data is available at present on 
the performance of the SES 100-B; either through computer 
simulation or actual craft tests. Since the computer 
program v;as installed at NPS, Monterey in October, 1972, 
this author was one of the first persons to use the program 
installed in the NPS 360/67* As a result, this author was 
like a growing lad; first learning to crawl, then walk, 
and finally to run as more and more knowledge was obtained 
regarding the computer program. 

B. .MULTIPLE CASES 

V/hen making multiple cases (runs) using one data deck, 
all control statements not specified for that particular 
case will carry the values of the specified data from the 
previous control statement over into the next case. 

C. LIMITATIONS 

The limitations that were determined [Ref. prior to 
use of the computer program were 

1 . Wave Components 

The maximum number of wave components that can be 
used is ten (10). 
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2. Wind 



No provision is made in the computer program for 
Inputing wind direction and wind speed. 

3 . Speed 

The minimum speed that can be used is the hump speed 
( 16.32 kts) vjhich is \^^here the craft goes into the ”on 
bubble” mode. 

4. V/ater Contact with Top of Plenum 

No provision is made in the computer program for 
calculating the effect of water striking the top of the 
plenum. The program does print the time. Immersion depth, 
and location along sldev/all where this v;ater contact does 
occur. 

5 . Draft 

The range of drafts must be selected so that the top 
of the plenum does not sink below the water surface. For 
the SES 100-B, a maximum draft of 6.17 feet is allovied. 

6. Fan Flow Rates and Fan Power 

Fan flow rates and the associated fan power required 
to produce the fan flow rates are considered to be instan- 
taneous values. No delay times are used in the computer 
program when fan power and fan flow rates are changed in 
value . 

D. INITIAL CONDITIONS 

The first data that must be obtained is the initial 
conditions for various speeds under calm water conditions. 



No waves are used so control statement 01102 (see Appendix 
B) is omitted. In addition the speed Is a constant input 
and the thrust Is allowed to vary by placing a one (1.0) In 
columns 26-35 of control statement 00105 (see Appendix B). 

To Insure that calm water steady-state conditions are 
reached, the problem time should run to at least 40 seconds. 
The initial steady- state variables that are required as 
output data are draft, pitch angle, plenum pressure, and 
thrust. This data can be obtained by placing a one (1) in 
block 15, 35, and 55 on the IBM card following control 
statement 00102 (see Appendix B). This data could be made 
part of the variables printed under the vertical and lateral 
plane summaries of control statement 00104 (see Appendix B). 
The initial condition values obtained are shown in Table I. 
The ratio of problem time to computer time was found to be 
1:0.8 for establishing initial conditions. 

To check to see if initial conditions are correct, a 
zero (0.0) is placed in blank 26-35, or leave blank, of 
control statement 00105. A blank is considered a zero in 
the computer problem. This establishes the value of thrust 
as initial condition and allows the speed to vary. The 
values of the initial condition variables should not change 
with time if the initial conditions are selected correctly. 

E. TOLERAIsICE VALUES OP INTEGRATORS 

The next item that must be established is the tolerance 
value to be assigned to each integrator. The recommended 
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TABLE I 



INITIAL CONDITIONS 

(AS DETERMINED UNDER CALM WATER CONDITION) 



SPEED 

(KTS) 


U 

(PT/SEC) 


DRAFT 

(FT) (IN) 


THETA 

(DEG) 


THRUST ONE 
(LBS) 


16.5 


27.87 


1.^420 


17. Oi] 


0.09 


8170.8 


20. 0 


33.78 


1.^06 


16 . 87 ' 


0.12 


6490.6 


25.0 


ii 2.22 


1.381 


16.57 


0.15 


5314.9 


30.0 


50.67 


1.351 


16.21 


0.18 


4882.7 


32.5 


5^.89 


1.33^ 


16.01 


0.19 


4832.4 


33.0 


55.73 


1.330 


15.96 


0.20 


4832.3 


35.0 


59.11 


1.315 


15.78 


0.21 


4860.2 


^ 0.0 


67.56 


1.272 


15.26 


0.23 


5091.4 


^ 5.0 


76.01 


1.223 


1 ^ 1.68 


0.26 


5490.2 


50.0 


8 ^.i <5 


1.166 


13.99 


0.28 


6002.1 


55.0 


92.90 


1.103 


13.2^4 


0.31 


'6589.3 


60. 0 


101.33 


. 1.033 


12 .i 40 


0.33 


7220.5 


65.0 


109.81 


0.956 


ll.i 47 


0.35 


7872.8 


70. 0 


118.22 


0 . 87^4 


10 .i ]9 


0.36 


8521.1 


75.0 


126.70 


0.786 


9.^32 


0.38 


9147.1 


80. 0 


135.11 


0.69^ 


8.33 


0.39 


9731.7 


85.0 


IM3. 61 


0.598 


7.176 


0.39 


10254.3 


90. 0 


152.0 


0 .i ^98 


5.98 


0 . il 0 


10697.3 


Note : 


Plenum pressure of 


92.80 PSP 


is used 


as initial 




condition 


for all 


speeds . 
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—5 

procedure [Ref. 4] was to set all values to the order of 10 
-9 

to 10 . Then half the value of one tolerance, rerun the 

program, and compare output data values. If the results do 
not change drastically, then double the tolerance again. As 
changes are noted in the integrator outputs, hold the toler- 
ances for those integrators fixed and vary other tolerances 
until the largest tolerances possible are obtained. 

The Inputs and output of each integrator are shovm in 
Fig. 1. Since simulation runs were made into or vjlth vjaves 
on a steady course, it was found that integrators number 
three (3) and ten (10) caused the greatest change in output 
data for a factor of ten (10) change in their tolerance. 

The selected values of tolerances for each Integrator 
that were used in all computer simulation runs are listed 
in Table II. 



TABLE II 

TOLERANCE VALUES FOR INTEGRATORS 



INTEGRATOR NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



TOLERANCE VALUE 

.0001 

.0002 

.000001 

.0000001 

.0001 

.000001 

.000000001 

.000000001 

.0001 

.0001 
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RKM = Runge-Kutta-Merson Method 
SR = Simpson's Rule 



Figure 1., Inputs and Outputs of Integrators 
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The values of the tolerances listed in Table II gave 
output data that had less than one (1) percent difference 
in values when the tolerance values were decreased by a 
factor of ten (10). If Integrator number three (3) or ten 
(10) v;as increased in value by a factor of ten (10), the 
output data would change by as much as six (6) percent. 

P. INITIAL PARAMETERS SELECTED 

The desired range of parameters initially selected to 
obtain the data requested by paragraph 2a [Ref. 2], are 
listed below. 

1. Speed 

Speed ranges of 30 knots to 90 knots were initially 
selected in 10 knot increments to be Investigated. 

2 . Wave Direction 

Waves directly ahead and astern were selected to be 
investigated. 

3 . Wave Length 

Wave lengths from one-half lamda (1/2X) to five 
lamda (5X) were initially selected to be investigated. X 
is equal to the length of the SES 100-B which is 72.0 feet. 

^ . Wave Amplitude 

Wave amplitudes from 1 foot to 6 feet were initially 
selected to be investigated. 

G. CHANGING WEIGHT 

Symmetry of weight distribution is assumed between port 
and starboard side of the SES 100-B when using the computer 
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program. A nodal method is used to specify the v;eight and 
center of gravity of each node. Control statement 00202 is 
used with the nodal method. Appendix C lists the various 
input values for the nodal distribution of weight and center 
of gravity. 

If additional vjeight is to be added to the craft, the 
value of the weight to be added should be halved since 
symjnetry is used. The value used as input will only repre- 
sent that amount added to the starboard side. The center 
of gravity specified will be the same distance to starboard 
as to port. 

Control statement 00201 uses the overall v;eight, center 
of gravity, and mass movement of inertia about x-axis, y- 
axis, z-axis, and xz-axis . 

Even though control statement 00201 uses only one IBM 
card, as compared to the numerous cards used with control 
statement 00202 (see Appendix C), it was found from numerous 
computer runs that the computer time decreased by approxi- 
mately 10 percent when control statement 00202 was used. 

H . PROBLEMS ENCOUNTERED 

1 . Minimum Speed 

The first problem encountered in using the computer 
program for the SES 100-B v/as that the longitudinal velocity 
(U) did not reach a steady-state value when waves v;ere used 
at a speed of 30 knots. Varying wave lengths from 36 
feet ( ) to 360 feet (51) and wave amplitudes from 
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0.01 feet to 1.0 foot at a speed of 30 knots caused the 
longitudinal velocity to decrease exponentially as l-e^^’^® 
until hump speed was reached. This result showed that the 
longitudinal acceleration did not go to zero in order for 
the longitudinal velocity (U) to reach a steady-state 
condition. 

An investigation into the various forces in the 

x-direction was conducted. Table III is a list of the calm 

water force components in the x-direction for various 

speeds. Force component PXPVJAV decays exponentially as 
-1 56 

e with an increase in speed. The other force compo- 

nents (PXBS, PXSS, PXSV/, PXRUD, PXAED) increase as with 
an Increase in speed. FXP and Thrust are equal and opposite 
in value and represent the sum of all the forces in the 
x-direction. A plot of FXPWAV and the other forces in the 
x-direction versus speed (Pig. 2) shows that the cross-over 
occurs at approximately 3^*0 knots. 

An investigation into the minimum critical speed 
(Uoc) was conducted (see Appendix A) v\rlth the result that 
Uoc was approximately equal to 31.0 knots. This result 
agrees fairly well with the cross-over of the e” and 

forces. Based on these results, the minimum speed used to 
obtain the output data was established at 35 knots. 

2. Maximum Wave Amplitude 

A 2 foot amplitude (height) wave is the same as a 
^ foot high wave measured peak to trough. Wave amplitudes 
of 6 and 8 feet were tried. In this particular case, the 
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TABLE III 



BREAKDOV/N OF FORCES 
(All Values in lbs) 



SPEED (KTS) 


FXBS 


FXSS 


FXSW 


FXRUD 


FXAED 


FXPWAY 


16.5 


172.4 


104.6 


795.6 


117.2 


116.3 


15036. 


20.0 


240.5 


148.3 


1130.2 


166.1 


170.8 


11125. 


o 

• 

C\J 


351.6 


221.6 


1694.1 


250.7 


266.9 


7845 . 


30.0 


476.4 


305.4 


2349.7 


352.0 


384.3 


5897. 


3 ^ 4. 0 


583.2 


378.9 


2933.3 


445.1 


493.7 


4847 . 


3^.5 


596.8 


388.5 


3009.5 


457.4 


508.3 


4738. 


35.0 


610.5 


398.0 


3086.3 


470.0 


523.1 


4632 . 


40.0 


749.0 


497.0 


3890.5 


604.3 


683.3 


3758. 


50.0 


1018.5 


702.1 


5644.1 


921.6 


1067.7 


2650. 


60.0 


1238.1 


891.5 


7478.9 


1302.9 


1537.7 


1991. 


70.0 


1357;4 


1028.8 


9250.4 


1747.5 


2093.5 


1564 . 


80.0 


1321.3 


1070.8 


10812.0 


2254.8 


2735.2 


1269. 


90.0 


1062.6 


968.0 


12023.0 


2823.6 


3462.1 


1055. 



0 

0 

0 

5 

6 

0 

5 

6 

3 

8 

n 

0 

3 
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Figure 2, Force Versus Speed for Calm Water Condition 



computer program v;ould stop In approximately 7 to 11 seconds 
of problem time due to the size of the incremental Integrator 
step size becoming smaller than 1 x 10*”^ seconds. Varying 
the tolerance values assigned to each integrator would not 
eliminate this condition. The result is that the Runge- 
Kutta-Merson method of integration is restricted to a certain 
range of V7ave lengths and wave amplitudes in order that the 
error generated, when compared with tolerance values, be such 
that the Incremental step size of integration would be halved, 
but not halved until the incremental step size v;as less than 
1 X lO"^ seconds. Some limits of speed, wave length, and 
wave amplitude are given in Table IV, for both ahead and 
astern waves, where the computer program stopped due to the 
incremental integrator step size becom.lng less than 
1 X 10”^ seconds. In every case, this condition occurred as 
a result of integrator number 3. 

3. Water Contact with Plenum 

The computer progarm will print out the time, immer- 

- N 

slon depth, and location along the sldev;all whenever the 
amplitude of the wave would strike the top of the plenum. 

No provision is made in the computer program to account for 
the additional forces created by the wave striking the top 
of the plenum. It v/as assum.ed that Immersion depths of only 
several inches above the top of the plenum v 7 ould have very 
little effect on the performance of the craft. Since the 
top of the plenum is 6.17 feet above the keel, immersion 
depths equal to or less than 6.6 feet were considered as 
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having negligible effect. Data obtained where the immersion 
depth was greater than 6.6 feet was considered invalid. 

Some limits of speed, wave length, and wave amplitude are 
given in Table IV for both ahead and astern vjaves. 

H . Negative Fan Horsepower 

Total fan horsepower is calculated based on the 
amount of fan inflovj rate into the bow seal, plenum, and 
stern seal. Under certain wave conditions, a negative fan 
inflow^ rate was indicated for either the bow seal, plenum, 
or stern seal. Normally the total fan inflow rate would 
range from 3,000 cubic feet per second to 10,000 cubic feet 
per second. Several different computer output runs indi- 
cated that the fan inflow rates to the bow seal v;ould vary 
between 1,000 cubic feet per second to -100 cubic feet per 
second. This small amount of negative fan inflow rate, when 
compared to the large net positive value of total fan inflow 
rate, v;as considered not sufficient to alter the output data. 
When the net total fan inflow rate was negative, which 
normally requires a negative fan inflow to all three loca- 
tions, a negative total fan horsepower resulted. Since the 
fans can not be put in reverse (pumping air out), this 
condition was considered severe enough to cause the output 
data to be invalid. Some limits of speed, wave length, and 
wave amplitude, for both ahead and astern waves, are listed 
in Table IV where negative fan horsepower resulted. 



35 



NON-USABLE PARAMETERS 
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5 . Problem Time Versus Computer Time 



The ratio of problem time to computer time for a 
1 foot amplitude wave v;as found to be approximately 1:^0. 
Increasing wave length and speed Increased this ratio very 
slightly (l:il2). This ratio holds good for both ahead and 
astern waves, with just a slight decrease in the ratio for 
astern waves (1:38). Increasing the v/ave amplitude from 
1 to 2 feet would change the ratio from 1:^0 to 1:60. 
Increasing the v/ave amplitude from 1 to 3 feet would change 
the ratio from 1:^0 to 1:70. 

6 . Reducing Computer Time 

It was found that using the calm water initial 
conditions would require approximately 1^0 to l60 seconds 
for the value of U to reach a steady-state condition. With 
an additional 10 to 20 seconds added for printing out steady- 
state values, a total of 2.0 to 3-0 hours of computer time 
was required for each combination of speed, wave direction, 
wave length, and wave amplitude. 

To reduce the computer time required, a run of ^0 
seconds problem time was used with the initial calm v/ater 
conditions as input values. The value of U obtained for 
each run was plotted between 0 and ^0 seconds and an esti- 
mate of steady-state speed (U) was made. This steady-state 
speed (U) could be estimated because the speed (U) would 
decrease exponentially with time from its calm water value 
as waves v;ere encountered. The new estimated steady-state 
speed (U) was then used in control statement 01201 as a new 
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initial speed, and the program was run for an additional 
30 seconds. This resulted in a reduction of 1.0 to 1.5 
hours of computer time for each parameter of speed, wave 
direction, wave length, and wave amplitude used. 

7 . Output Data Limitation 

The IBM- 360/67 computer at the V/.R. Church Computer 
Center has a built-in limitation as to the maximum number 
of pages of output data that can be obtained. A rule of 
thumb used v;as that a maximum of 100 time intervals of 
data could be printed out if all the print switches listed 
under control statement 00102 were ON (indicated by a 1.0). 
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IV. PRESENTATION OP DATA 



A. DRAG VERSUS SPEED 

A plot of drag versus speed under calm water conditions 
is shovm in Pig. 3. The results show that a minimum value 
(saddle point) in the curve occur’s at a speed between 32.5 
and 33*0 knots. Comparing Pig. 3 with the results reported 
in Ref. 5j close agreement does exist between the value of 
hump speed and the value of drag occurring at hump speed. 
Pigure 3 shows that the minumum value of drag occurs at a 
lower speed and that this particular value of drag is less 
than that reported in Ref. 5. 

Pigure 3 indicates that once hump speed is reached and 
the craft goes "on the bubble", it would accelerate from 
16.5 knots to approximately 68 knots under calm vjater 
conditions if the value of thrust was held constant. 

Reference 5 indicated that the craft would accelerate from 
hump speed to a speed of approximately ^0 knots if the 
value of thrust was held constant. 

B. STEADY-STATE SPEED 

Each particular wave condition resulted in a different 
steady-state speed being obtained. Table V lists the steady- 
state speeds that were obtained from various wave directions, 
wave lengths, and v.^ave amplitudes. Although the steady-state 
value of speed is only valid for the particular wave condi- 
tion listed, the values are useful in estimating other 
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Figure 3 « Dreg Versus Speed for Calm Vlater Condition 







<D 










































C 4-> ^ 
O CO 

•H 4-^ E*H 




C\J 


I — 1 




CM 


CO 


iH 


-=r 


CO 


t-- 


CO 


CM 


VO 




CM 


CM 


CO 


VO 






+i t/j w 


o 


o 


C\J 


o 


O 


o 


CM 


o 


o 


co 


o 


o 


I — 1 


o 


o 


o 


r — 1 


o 






cd 1 ^ 
>> 

4^ ^ 

O Clj CL) 
^ 0) CD 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 






r-H 4-^ 

pq CO CO 










































State 

KTS) 










































1 


•=d* 


o 


I — 1 


'O 


CM 


1 — 1 




CM 


m 


VO 


m 


o 


m 


CO 


o^ 


LTV 


o\ 


CO 






>> 

^ X3 


CO 


t>- 


CO 


CO 


1 — 1 


VO 


-=r 


cr\ 


o 


•=4* 


CO 


CM 


cr\ 


Lpy 


CM 


VO 


IT\ 


CO 






Cj (D 


o 




CO 


l>- 


VO 


CTN 


LTV 


CM 




CO 




Or\ 


CM 


CO 


CM 


c^ 


G\ 


CO 






(1) CD 
-P Qa 
CO CO 




CO 


Cvl 


C\J 


CO 


CO 


CO 


CO 


CO 


CM 








CO 


m 


m 




















































Eh 










































pH 








































CO 










































Q 










































W 


CD 








































W 


"d 


















































































CO 


4^ 






































> 




*H 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 




w 


CD iH 






































w 


Eh 


> a 


1 — 1 


1 — 1 


! — 1 


CM 


1 — 1 


1 — 1 


1 — 1 


CM 


CM 


CM 


fH 


1 — 1 


1 — 1 


CM 


1 — ! 


! — 1 


1 — 1 


CM 




< 


Cd S 






































CQ 


Eh 


:b: < 






































< 


CO 










































1 

>H 










































Q 


Eh 








































< 


Ph 








































W 










































Eh 










































CO 


x: 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 






4^ 










































CD bO 


VO 


C\J 


oo 


MD 


VO 


CM 


CO 


VO 


CM 


CO 


VO 


CM 


CO. 


VO 


VO 


CM 


CO 


CM 






> d 


CO 


tr^ 


o 


CO 


CO 


tr^ 


o 


CO 




o 


CO 


tr^ 


o 


CO 


(O 


ir^ 


o 








cd CD 

:s: ^ 






1— 1 








I — 1 






I — 1 






I — 1 








I — 1 








c 










































O 










































•H 






'O 






x^ 


x^ 


X3 




x^ 


x^ 


X3 


x^ 








x^ 


x^ 






4^ 


cd 


cd 


cd 


Cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 


cd 






O 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 






CD CD 






x: 




x: 




x: 








XI 












x: 








> Pi 
cd *H 
Q 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 
















































CO 










































Eh 




















































































I — 1 

cd 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 






•H ^ 
4^ CD 


m 


LTV 


m 


lf\ 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 






•H CD 
C P. 
M CO 


CO 


CO 


CO 


CO 


•=d* 


-=d* 


-=r 


-=d* 


-=4* 




LTV 


m 


ir\ 


in 


VO 


VO 


VO 


VO 



ill 



Initial Wave Wave Wave Steady-State Fluctuation 

Speed (KTS) Direction Length (FT) Amplitude (FT) Speed (KTS) Steady-State 

Speed (KTS) 



CM 




CM 




•=T 


m 


\o 


on 


in 


cr\ 


CM 


m 


CX) 




CT\ 




m 






O 


O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


CM 


o 


I — 1 


o 


o 


1 — 1 


1 — 1 


O 


O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 



cr\ 


cn 


o\ 


CXD 


MD 




CM 


in 




CM 


cn 


vn 






MD 


CXD 


CT\ 


CXD 


in 


H 




KO 


CT\ 


MD 


lf\ 


m 








m 


I — 1 


o 








cr\ 




-=r 


rH 




VD 


m 




CT\ 


CXD 


CT\ 


MD 


o\ 


o 


CM 




cn 




o 


o 


MD 


\o 


MD 




cn 


cn 


CM 


cn 




cn 










cn 


m 


in 


in 


in 



o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


1 — 1 


1 — 1 




1 — 1 


1 — 1 


I — 1 


CM 


1 — 1 


1 — 1 


CM 


1 — 1 


1 — 1 


1 — 1 


CM 


cn 


1 — 1 


1 — 1 


1 — 1 


1 — 1 



o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


\o 


CM 


MD 


CM 


VD 


CM 


MD 


MD 


CM 


MD 


MD 


CM 


CXD 


MD 


MD 


MD 


CM 


CXD 


cn 




cn 


1-- 


cn 




cn 


cn 




cn 


cn 




o 
1 — 1 


cn 


cn 


cn 




O 
1 — 1 















C 


C 








d 


C 


c 




d 


C 


C 








T5 




u 


u 












P 


P 


p 


u 


U 


U 


P 




cS 


cd 


cd 


cd 


CD 


<D 


CD 


CD 


<D 


<D 


CD 


CD 


CD 


(D 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


CD 


4^ 


-P 


4^ 


43 


43 


43 


43 


43 


43 


43 


43 


43 


43 


43 


43 




x: 


ja 




W 


W 


W 


W 


W 


W 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 


< 



o 


o 


o 


O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


O 


O 


in 


in 


o 


o 


in 


in 


in 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


O 


O 






<D^ 


CD^ 


cn 


cn 


cn 








in 


in 


in 


in 


in 


MD 


MD 


MD 


MD 



k2 



Initial Wave Wave Wave Steady-State Fluctuation 

Speed (KTS) Direction Length (FT) Amplitude (FT) Speed (KTS) Steady-State 

Speed (KTS) 



OJ 




j=r 


VO 




oo 


CM 


o 




o 


o 


o 


CM 


o 


1 — 1 


oo 


o 


o 


o 


o 


O 


o 


o 


o 



o 


VO 


I— 1 


C-- 


oo 


CM ^ 


oo 


CM 


OO 


1 — 1 


vo 


CM 




o 


t — 1 


CM 


o 


OO 


oo 


oo 


cta l>- 


o 


Lf\ 


LPi 


VO 


VO 


VO 


VO 


t-- 





oooooooo 



o 


o 


o 


o 


o 


o 


o 


o 


VO 


VO 


VO 


CM 


oo 


VO 


CM 


oo 


1—1 

CM 


oo 


oo 




o 
1 — 1 


CO 


t-- 


o 

1—1 



CD 

-P 

to 



C 

(D 

P> 

to 

< 



c 

Sh 

<L) 

to 



a 

(D 

to 

< 



a 

<L) 

to 



(D 

to 



C 

CD 

-P 

to 

< 



C 

<D 

P> 

to 

< 



o 


o 


o 


o 


o 


o 


o 


o 


o 


o 




LT\ 


m 


o 


o 


o 


VO 


VO 










CJ^ 


CT\ 



^3 



steady-state speeds under vjave conditions that are not 
listed. This will reduce computer usage time if the proce- 
dure listed in Paragraph III H-6 is followed. 

C. BOW AND SIDEWALL LEAKAGE 

Under certain vmve conditions, air leakage would occur 
from the bow and/or sidewalls. Of course this is in addi- 
tion to the stern leakage that occurs under all sea 
conditions. To better understand when leakage from the bow 
and/or sidewall did occur. Table VI lists the steady-state 
craft speed, wave direction, viave length, and wave amplitude 
when this leakage occurred. 

D. PHASE RELATIONS OF VARIABLES 

To help understand hov; the simulated motions of the 
craft would respond to varying wave conditions, a phase 
comparison of selected variables with respect to wave motion 
was conducted. It was found that changing the .speed of the 
craft did not alter the phase relation betv;een the selected 
variables and the particular wave used under both ahead and 
astern wave directions. In addition, changing the v/ave 
amplitude did not alter the phase relation between the 
selected variables and the particular wave used under both 
ahead and astern wave directions. A representative initial 
condition speed of 60 knots was selected to show the phase 
relation between the different variables with respect to 
various v;ave lengths. 
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TABLE VI 

OCCURRENCE OF BOW AND/OR SIDEV/ALL LEAKAGE 
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TABLE VI (Continued) 
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The variables that have been selected for study, based 
on Ref. 2 are: Draft, Plenum Pressure, Bovj Acceleration, 

Center of Gravity Acceleration, Bow Seal and Sidcv/all 
Leakage Rates (where applicable),' Stern Seal Leakage Rate, 
Input Fan Plow Rates (bow, plenum, stern), and Total Pan 
Horsepov/er . 

1. Ahead V/aves 

Wave lengths of 3^ feet (1/2 X), 72 feet (X), and 
108 feet (3/2 X) are shovm in Figs. ^ through 2k. 

Figures 5, and 6 shov/ that the response of the 
center of gravity acceleration leads the wave for a wavc“ 
length of 36 feet and lags the v/ave for v;ave lengths of 72 
and 108 feet. The smallest magnitude of center of gravity 
acceler-ation occurs for the 72 foot V'lavelength . Figures 
7, 8, and 9 show that the bov/ acceleration is in phase v;ith 
the center of gravity acceleration and its magnitude is 
larger. This is to be expected for a rigid type ship. 

The average draft response is shovm in Figs. 10, 11, 
and 12. No distinct phase relationship exists since the 
values of draft are averaged along the length of the side- 
v;all. Plenum pressure is the primary determinant in estab- 
lishing average draft values. The higher the plenum pressure, 
the smaller the average draft value that will occur. The 
72 foot wave length has the smallest draft with the corre- 
sponding largest plenum pressure (see Pig. l4). 

Figures 13, 1^, and 15 show the response of plenum 
pressure to different wave lengths. Pressure in the plenum 
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is primarily based on the plenum volume and the difference 
between the air flow in and out of the plenum. As the 
volume is decreased (Increase in draft), the plenum pressure 
will increase. Comparing plenum pressure with average 
draft, you can see that the two are in phase. In addition, 
the pressures into the bow and stern seals are in phase 
with the plenum pressui’e. 

Figures l6, 17, and l8 show the stern seal leakage 
that occurs under different wave lengths. The value of 
stern seal leakage is based upon a calm v;ater baseline 
leakage rate of 6^117 cubic feet per second at a plenum 
pressure of 92.8 pounds per square foot. Since it is not 
known precisely V'jhere all the leakage does occur in the 
craft, it is assumed that all of the baseline leakage does 
occur out the stern seal. The values of stern seal leakage, 
under wave conditions, are based upon the principle of 
Venturi flov;. For the same area of stern seal gap, the 
amount of stern seal leakage will vary proportionally with 
the plenum pressure. As a result, the stern seal leakage 
rate is in phase v/ith the plenum pressure. No bow or 
sidewall leakage occurred at the wave lengths shov;n. 

The input fan flow rates to the bow seal, stern 
seal, and plenum are shovm in Figs. 19, 20, and 21. The 
fan flov; rates are based on a function of pressure. Since 
the fan characteristics are nonlinear, a series of linear 
curve-fits for the pressure versus. input fan flow rate is 
used. The slope of the pressure versus input fan flow rate 
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changes at different locations on the plot, with the result 
that the fan flow rate causes a more rapid increase in fan 
flow rates at lower pressures. Basically the fan flow Inputs 
are out of phase with the pressures. In Pig. 19 and 21, 
negative fan flov; Inputs to the bow seal are not shovm. 

The maximum negative value that did occur v;as approximately 
-I 50 pounds per square foot. 

Figures 22, 23, and 24 show the total fan power 
required to produce the necessary input of air flow. Of 
interest here is to note the double frequency of oscillation 
that occurs. The total fan power is based upon a summation 
of fan power required for the bow seal, plenum, and stern 
seal. In turn each individual fan power is based upon the 
pressure times the fan flow input. Since pressures and fan 
flow Inputs are out of phase, the maximum power requirements 
will occur when the pressures are the largest and also v.^hen 
the fan flov; inputs are the largest. This results in the 
double frequency of oscillation of the total fan power 
requirement. 

2 . Astern Waves 

Wave lengths of 72 feet (X) and 108 feet (3/2 X) 
are shown in Figs. 25 through 40. Wave lengths of 36 feet 
are not shown for astern vmves since they follov; approximately 
the same sinusoidal phase relation established for the 36 
foot wavelength of ahead waves. 

Figures 25 and 26 show the relationship of the 
center of gravity acceleration to the viave. The maximum 
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acceleration lags the v;ave by a small amount. Figures 27 
and 28 show that the bow acceleration is in phase with the 
center of gravity acceleration, as is to be expected for 
a rigid type ship. 

The average draft response is shown in Fig. 29 and 
30 . No phase comparison between draft and wave can be made 
since this is an average value. The value of draft is 
primarily dependent on the value of plenum pressure. 

Figures 31 and 32 show the response of plenum 
pressure to different vjave lengths. Basically the plenum 
pressure follov.'s the average draft response, but Fig. 32 
shows hov/ the plenum pressure can deviate from a sinusoidal 
motion. Doth bov/ seal and stern seal pressures are in 
phase with the plenum pressure. 

Figures 33 and 3^ show the stern seal leakage for 
different wave lengths. The stern seal leakage rates are 
in phase with the plenum pressure response as a result of 
the same reasoning used for ahead waves. 

Figures 35 and, 36 shov; the bow and sidev;all leakage 
rates. No sidewall leakage occurs at the 108 foot wave length 
(Fig. 36 ). Leakage will only occur when a gap is created 
between the bow seal and water surface and/or the sidev/all 
and water surface. The amount of leakage rate is dependent 
upon the area of the gap created and the pressure of the 
plenum. The principle of Venturi flow applies to the leak- 
age rate. V/hen sidevjall leakage did occur (Fig. 35), this 
leakage occurred at the forv/ard part of the sidewall. 
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The Input fan flov; rates are shown In Figs. 37 and 
38 . As v;lth ahead waves, the input fan flov; rates are out 
of phase with the pressures. 

Figures 39 and ^10 show the total fan power required. 
Again notice the double frequency of oscillation that occurs. 
This occurs because of the same reasoning used for explaining 
the ahead wave case. 

E. AVERAGE AND FLUCTUATION VALUES OF VARIABLES 

Average and fluctuation values of selected variables 
for various speeds, wave directions, v;ave lengths, and v;ave 
amplitudes are shown in Figs. ^1 through 80 . The variables 
that are shown In the different figures are: Pitch, Center 

of Gravity Acceleration, Bow Acceleration, Draft, Bov; and 
Stern Seal Pressures, Plenum Pressures, Stern Seal Leakage 
Rate, Bow Seal Input Pan Plow Rate, Plenum Input Pan Plow 
Rate, Stern Seal Input Fan Flov; Rate, and Total Pan Power. 

The fluctuations, and not the average values, of Bov; Acceler- 
ation and Center of Gravity Acceleration are shov;n since 
their average values are approximately zero. 

The average values were calculated by adding the values 
of a number of data points (20 to 50) along a full wave 
encounter cycle and dividing by the number of data points. 
Either one or more full v^ave cycles vrere used depending on 
the number of data points available in each wave cycle. 

This summation of values v;as necessary because of the 
non-slnusoidal nature of certain variables. 
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The values of fluctuations used v/ere based on taking 
the difference between the maximum and minimum values of 
the variables. 

The dots shown In Pigs. ^1 through 80 represent the 
data points. No Interpolation was used between the data 
points . 

The values of vjave length and wave amplitude are denoted 
on the figures at the right hand side of the last data point 
used. A 36 foot v/ave length and 1 foot wave amplitude is 
denoted as 36 x 1 . 

Figures ^11 through 60 are for ahead waves (directly on 
the bow) and Figs. 6 l through 80 are for astern v/aves 
(directly astern). 

A discussion of the results obtained in Figs. ^1 
through 80 is contained in later sections. 
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Figure 4 • Response of Center of Gravity Acceleration to Ahead Wave 
Initial Conditions: Speed 60 Knots, Wave Length 36 Feet, 'nave Amplitude 1 Foot 
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Figure 5 ♦ Response oi Center ojr Gravity Acceleration to Ahead V^ave 
Initial Conditions: Speed 60 Knots^ Vfeve Length 72 Feet, VJave Amplitude 1 Foot 
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Figure 6 , F.esponse oi Center of Gravity Acceleration to Ahead vi'ave 
Initial Conditions: Speed 60 Knots j VJave Length 103 Feet, Vfave Amplitude 1 Foot 
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Figurs 7* Response of Bo'ix Acceleration to Ahead Vfeve 
Initial Condit5.cns; Speed 60 Knots, Have Length 36 Fest, Vlave Amplitude 1 Foot 
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Figure ® « Response of 3ovr Acceleration to Ahead V/ave 
Initial Conditions; Speed 60 Knots^ Wave Length 72 Feet, Vfeve Amplitude i Foot 
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Figuro 9 • Response of Bow Acceleration to Ahead V/ave 
Initial Conditions; Speed 60 Knots, Vfave Length 103 Feet, Wave Amplitude 1 Foot 
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Figure 10, Draft Rseponse to Ahead Wave 
Initial Conditions: Speed 60 Knots, Wave Length 36 Feet^ Wave Anplitude 1 Foot 

D-raft 1.033 Feet 
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Figure X3« Response oS Plenisa Pressure to Ahead Wave 
Initial Conditions: Speed 60 Knots, V/ave Length 36 Feet, V7avs Amplitude 1 Foot, Plenum Pressure 92. S PSF 
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Pleniaa Ptcsstito 92 «S PS? 
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Figure 15 « Response of Plenxsn Pressure to Ahea.d Uave 
Initial Conditions; Speed 60 Knots, Vfave Length 103 Feet, Vfave Amplitude 1 Foot 

Plenum. Pressure 92,8 PSF 
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Fisvjro 16 a Stem Seal Lsalcage Rate for Ahead Wave 
Initial Conditions; Speed 60 Knots.. Wavo Length 36 Feet, Wave Amplitude 1 Foot 

Stem Seal Leakage Eate 6417 Cubic Feot/Socond 
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Figure IS , Stem Seal Leakage Fate for Ahead Wave 
Initial Conditions; Speed 60 Knots, Vfeve Length 103 Feet, Vfeve /mplitude 1 Foot 

Stem Seal Leakage Rate 6417 Cubic Feet/Second 




(03S/ii no) MOTi VKi ^Tidui (jj) *0 Jo 'o ©pnqfidnnr 



68 



Figure 19 • Input Fan Flow Rates for Ahead Wave 
In5.tial Conditions: Speed 60 Knots, Wave Length 36 Feet, Wave Amplitude 1 Foot 
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Figure 20, Input Fan Flew Rate a for Ahead Viave 
Initial Conditions; Speed 60 Knots, Wave Lenfjth 12 Feet, Waw Amplitude 1 Foot 
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Figura 21, Input Fan Flo^r Rates for Ahead Vjavc 
Initial Conditions; Speed 60 Knots, Wave Length lOS Feet, Vi'ave Amplitude 1 Foot 
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Figure 22, Total Fan Po^ver Response to Ahead VJave 
Initial Conditions: Speed 60 Knots, Wave Length 36 Feet, VJave Amplitude 1 Foot 
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Figure 23 » Total Fan Power Response to Ahead Wave 
Initial Conditions: Speed 60 Knots, VJave Length 72 Feet, V/ave Amplitude 1 Foo'J 
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Figure 24« Total Fan Povrsr Resj>onGO to Ahead Vlavc 
Initial Conditions: Speed 60 Knots, Vfevo Length ICS Feet, V/ave Amplitude 1 Foot 
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Figure 25« Response of Center of Gravity Acceleration to Astern VJavo 
Initial Conditions; Speed 60 Knots, V/ave Length 72 Feet, Wave Amplitude 1 Foot 
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Figviro 26* Response of Center of Gravitj Acceleration to Astern V<ave 
Initial Conditions: Speed 60 Knots^ Wave Length lOS Feet;, Wave Amplitude 1 Foot 
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Figure 27 » Response of Bow Acceleration to Astern Wave 
Initial Conditions: Speed 60 Knots, Vfave Leng'bh 72 Feet, Wave Amplitude 1 Foot 
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Figure 29 e Draft Response to Astern VJave 
Initial. Conditions: Speed 60 Knots^ VJave Length 72 Feetj, VJave 

Amplitude 1 Foot, Draft 1,035 Feet 
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Figurs 30, 35raft Response to Astern V/ave 
Initial Conditions: Speed 60 Itnots, Wave Length lOS Feet, V/ave 

Anplitude 1 Foot, Draft 1,033 Feet 
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Figure 34 e Stera Seal Leakage Rate for Astern Vfave 
Initial Conditions: Speed 60 Knots, vv'ave Length 103 Feet, Vfeve Amplitude 1 Foe 

Stem Seal Leakage Rate 6A3.7 Cubic Feet/Socond 
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Figure 35 • Eov; Seal and Sidevw’all Leakage Rate Tor Astern Wave 
Initial Conditions: Speed 60 Knots, Vlave Length 72 Feet, Vfeve Amplitude 1 Foo' 
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Figure 36 , Bo>r Seal Leai-3.ge Rate for Astern Wavs 
Initial Condition: Speed 60 Knots, Wavs Length 108 Feet, War^fe Amplitude 1 Foot 
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Figure 37 * 3i:put Fan Flow Pates for Astern \7avo 
Initial Conditions: Speed 60 Knots, Wave Length 72 Foet, Wave Amplitude 1 Foe 
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Figure 33 o Inpu't Fan 71c\f Rates for Astern Wave 
Initial Conditions; Speed 60 Knots, Viave Length 108 Feet, Wave Amplitude 1 Foot 
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Fxguro 39 * Total Fan Power Response to Astern V7ave 
Initial Conditions: Speed 60 Knots, V/ave Length 72 Feet, 'dave Amplitude 1 Foot 
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Figure 40 • Total Fan Power Response to Astern Wave 
Initial Conditions: Speed 60 Knots^ Vfeve Length 108 Feet, Wav® Amplitude 1 Foot 
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Figure Ul, Average Fitch Angle Versus Steadj^’-State 
Speed for Different Ahead Wave Cenditiens 
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ire lf2„ Fl’actuation of Pitch Angle Versus Steady-State 
Speed for D3-fferent Ah.ead Wave Conditior^ 
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Figure io. Fluctuation of Center of Gravity Acceleration Versus 
Steady-State Speed for D:lfferent Ahead V/avc Condj,tic>n3 
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,’xo Fluctuation of Eoir Acceleration Versus Stcady-Stat' 
Speed for Different Ahead V.'avs Conditions 
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Figure 45 «• Average Draft versus Steady-^tate Speed 
for Different Ahead Wave Conditions 
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Figure 46. Fluctuation of Draft 'Versus Steady-State 
Speed for Different Ahead V/ave Conditions 
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Figure 47* Average Bow Seal and Stem Sea?. Pressure Versus 
Steadj-^State Speed for Different Ahead Vfeve Conditions 
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Figure i>8, Flx\Cwuation of Bow Seal and Stem Seal Pressxxre Versus 
Steady-State Speed for Different Ahead V^a^-e Conditions 
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iguro 49* Average Plcnina Fressure Versus Stsad^*-State 
for Different Aliead Wave Conditieixs 
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50o Fluctuation of PD.enun Prescure Versus Steaciy-State 
Speed for Different Ahead V*ave Conditions 
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Figure 51 ♦ Average Stem Seal Leakage Rate Versus Steady-State 
Speed for Different Ahead VJave Conditions 
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Figure 5^* Fluctuat5.on of Stem Seal Leakage Rato Versus Steadjr-»State 
Speed for Different Ahead Wave Conditions 
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Figure 53 « Average Bov; Seal Input Fan Flov; Pate Versus Steadjr-Stato 
Speed for Different /'.head V7ave Conditions 
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Figure 54* Fluctuation of Bov; Seal Input Fan Flow Rato Versus 
Steady-State Speed for DiiTerent Ahead Viavo Conditions 
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Figure 55 « Average Plemm Biput Fan Flo'-r Rato Versus Steaeij'-Stato 
Speed for Different Ahead Wave Conditior^s 



CU FT/SEC 




§ 

- 4 ^ 






O 

O 

o 

C\i 



8 

o 



105 



Figure 56# Fluctuation of Plenum Input Fan Flow Rate Versus Steady-Stato 
Speed for Different Ahead V7ave Conditions 
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Figure 57, Average Stem Seal Input Fan Flo’.'r Rate Versus Steady-State 
Speed for Different Ahead V7ave Conditions 
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Figure 58. Fluctuation of Stem Seal Input Fan Flow Rate Versus 
Steady-State Speed for Different Ahead 77ave Conditions 
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Figuro 59» Average of Total Fan Po'wer Versus Steady-State 
Speed for Different Ahead Vfeve Conditions 
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Figure 60* Flucti:iation of Total Fan Power versuc Steady-State 
Speed for Different Ahead Vfeve Conditions 
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Figure 61 Average Pitch Angle Versus Steady-State Speed 
for Different Astern VJavo Conditions 



DEG 



rH 




CO 



lA 



o 

u\ 



tA 



o 



lA 

H 




O 



111 



Figure 62« Fluctuation of Fitch Angie Versus Steady-State 
Speed for Different Astern Vfeve Conditions 
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Figure 63. Fluctuation of Center of Gravity Acceleration Versus Steady-State 
Speed for Different Astern VJave Conditions 
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Figure 64. Fluctuation of Bow Acceleration Versus Steady-Sta,tc 
Speed for Different Astern Wave Conditions 
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FigTire 65. Average Draft Versiis Steady-Stats Speed for Different 

Astern Wave ConditLoriS 
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'igure 66, Fluctuaticn of Dra.ft Versus Stsady-Ste- 
Speed for Different Astern V/ave Conditions 
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Arerage Bow Seal and Stem Seal Pressure Versus Steady. 
Speed for DiAferent Astern Wave Conditions 
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Figure 68, Flxictue-tion of Bow Ses.1 and Stem Seal Press'ore Versus 
Steady-State Speed for Different Astern Wa\-e Conditions 
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Figure 69^ Averags Plenixn Pressure Versus Steady-State Spe: 
for Different Astern VJave Conditions 
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for Daffcrent Astern Wave Condition; 
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Figure 71 » Average Stem Seal Leakage Rato versus Steady-Sta- 
Speed for DilTcrent i\stem Vfave Conditions 
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Figure 72* Fluctuation of Stem Seal Leatege Rats Versus Steady-State 
Speed for Different Astern V/ave Conditions 
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73, Average Bow Seal Input Fan Flow Hate Versus Steady— State 
Speed for Different Astern Wave Conditions 
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Figure 74* Fiiictuation of Bow Seal Input Fan Flow Rate Versus Steady-State 
Speed for Different Astern VJave Conditions 
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Figure O, Average Plenum Input Fan Flow Rato 'versus Steady-State 
Speed for Different Astern Wave Conditions 
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Figure ?6» Fluctuation of Plenua Input Fan Flow Rate Versus Steadj-Sta' 
Speed for Different Astern Wave Conditions 
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Figure 77 ® Average Stem Seal Input Fan Fio'-v* Rate Versus Stcady-Stat' 
Speed for Diffei'ent Astern Wave Cenditions 
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Figure Fluctuation of Stem Seal Input Fan Flov; Rate Versus 
Steady-»State Speed for Different Astern V/ave Conditions 
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Figuro 79* Avers-ge of Total Fan Power Versus Steady-State Speed 
for Different Astern Wave Conditior^ 
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for Different' Astern v/ave Condi' 



V. DISCUSSION 



A. CAPTURE AIR BUBBLE CONCEPT 

The Capture Air Bubble (CAB) ship is designed to run 
along the viater surface on a cushion (bubble) of air. The 
concept is to "capture" the air supplied by fans in the 
plenum by using rigid sidev;alls and flexible bow and stern 
seals. The air that is captured in the plenum area will 
reduce the draft of the craft and in turn will reduce the 
amount of drag. This allows a higher speed at a given 
amount of thrust for the CAB ship than a conventional 
monohull craft of the same gross weight. 

In the computer simulation of the lOO-B SES, an initial 
base leakage of 6^17 cubic feet per second was assumed to 
occur out the stern seal under calm vrater conditions at a 
plenum pressure of 92.8 PSF. The plenum volume, under calm 
water conditions, varied betv/een 9,200 and 10,200 cubic feet. 
The variation in plenum volume is due to the speed of the 
craft. At slower speeds the draft was greater causing a 
decrease in plenum volume. 

No matter what the speed of the craft is, the calm v-rater 
results indicated that the volume of the plenum was replenished 
approximately every 1.5 to 2.0 seconds. A leakage of 6^17 
cubic feet per second may not be considered as extremely 
excessive, but the fact that the plenum is replenished 
every 1.5 to 2.0 seconds appears to be in contradiction to 
the theory of the "captured" air bubble. 
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Flgui^es 51 and 71 shov; that Increasing speed and increa- 
sing the wave amplitude v;ill reduce the average rate of 
stern seal leakage. The average plenum pi*essures and the 
bow and stern seal pressures also decrease with an increase 
in speed and v;ave amplitude as shown in Figs. My, ^9, 67, 
and 69 . This agrees v;ith the computer program since leakage 
rate is dependent upon the area of the gap created and the 
associated pressures. However, the average input of air 
into the bov; seal, stern seal, and plenum also decreases vjlth 
an increase in speed and wave amplitude as shown in Figs. 

53, 55, 57, 73, 75, and 77. This is contrary to the computer 
program where a decrease in pressure causes an increase 
in input air by the fans and vice versa. The average total 
fan povjer required also decreases v;ith an increase in speed 
and wave amplitude as shown in Figs. 59 and 79. 

The computer program assumes an initial gap at the stern 
seal of 25.52 square feet. An additional gap is calculated 
due to wave conditions. This additional gap v;lll occur only 
when the wave surface is approximately below the level of 
the keel. The total gap at the stern is then based on the 
initial gap and any gap due to the wave. It was found that 
the area of the gap due to the wave had little or no effect 
on the total area of the gap at the stern. This was pri- 
marily due to the small amplitude waves used in the study. 

From the results obtained of the input of air by the 
fans and the total fan power required, it shows that the 
more severe wave conditions require less input of air and 
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horsepower. The volume of the plenum vms still being 
replenished by air about every 2.0 seconds. 

B. WAVE CONDITIONS 

V/ave lengths of 72 feet (X), both for ahead and astern 
waves, had the smallest effect on the performance of the 
100-B SES. The amount of decrease in speed from the initial 
speed input to the steady-state speed v/as the smallest for 
a 72 foot v\rave length. In addition, the center of gravity 
acceleration and bow acceleration were less for the 72 
foot wave length than either the 36 foot or IO 8 foot wave 
length. Generally the average stern seal leakage rate, 
total fan power, bow' seal input of air, plenum input of 
air, stern seal input of air, and bow seal, stern seal, 
and plenum pressures vjere higher. 

Although it is not knovm why this result did occur, it 
appears that a wave length approximately equal to the 
length of the craft gives the best performance of the 
craft in regard to minimum motion. 
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VI. CONCLUSIONS 



A. PERFORMANCE 

The performance of the SES 100-B, as determined by 
computer simulation, did not vary appreciably between ahead 
(directly on the bow) waves and astern (directly astern) 
waves under the same wave conditions. Steady-state speeds 
were larger, average plenum pressures were slightly higher, 
the fluctuation of pitch angle was larger, and the average 
draft v;as slightly less for astern v^aves as compared to 
ahead v/aves. The other performance variables shown in 
Pigs. ^1 through 80 were of similar values. 

Under similar wave conditions, the computer program 
showed that the craft could proceed through larger amplitude 
astern waves than ahead waves. Table IV showed, for 
example, that a wave length of 36 feet and wave amplitude 
of 2 feet from astern could be' run at a speed of 60 knots 
(vice only 50 knots if the wave was from ahead) before a 
negative fan power developed. The craft could run at a 
speed of knots in an astern wave of wave length 36 feet 
and wave amplitude of ^ feet before an immersion depth 
greater than 6.6 feet occurred. At the same speed, a 
negative fan povrer condition developed for an ahead wave 
of • wave length 36 feet and v;ave amplitude of 3 feet. V/ave 
lengths and wave amplitudes in this range are associated 
with a state 2 to state 3 sea condition. Therefore, it 
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appears that the CAB concept, and particularly the SES 100-B, 
will be limited to calm or moderate sea conditions. 

When bow and sidewall leakage did occur, generally the 
amount of leakage was larger at the slov/er speeds for a 
given wave condition. Table VI shows when bov; and sidevmll 
leakage did occur. Increasing the wave length resulted in 
larger leakage rates out the bow and sidewall. Under a 
given speed and wave condition, the magnitude of bov; and/or 
sidewall leakage was larger for an astern wave than an 
ahead wave. 

As speed v;as increased, in general the average values 
of pitch and draft increased v.flth a corresponding decrease 
in their fluctuation values. The fluctuations of the center 
of gravity acceleration and bov/ acceleration also increased. 

As speed was Increased, in general the average values of 
bow seal, stern seal, and plenum pressures decreased v;lth 
an increase in their fluctuation values. The average 
values of stern seal leakage also decreased, in general, 
with a corresponding increase in fluctuation as speed was 
Increased. The average Inputs of air by the fans to the 
bow seal, plenum, and stern seal decreased with a corres- 
ponding Increase in fluctuation as speed was increased. 

The same trend also held for the average and fluctuation 
values of total fan pov/er. 

B. RECOMMENDATIONS 

The follov/lng recommendations of additional areas of 
study are made on the basis of the results obtained during 
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this study. 

1 . Leakage Ratos 

The I'esults of the computer program shovved that the 
air is being replenished in the plenum every 1.5 to 2.0 
seconds. From the CAB concept, this appears to be an 
excessive amount of air input required. It is recommended 
that the results of air leakage be compared against those 
of the SES 100-B to see if this is actually the case. 

2 . Air Inputs by Fans 

The results of the computer program showed that an 
average decrease in air input by the fans and an average 
decrease in the amount of fan pov.'er occurred v/hen a more 
severe wave condition v.^as encountered. This does not 
appear feasible from the CAB concept. It is recommended 
that a study of the modeling of the air inputs by the fans 
and leakage rates be conducted to resolve this apparent 
discrepancy. 

3 . Stern Seal Leakage 

The amount of stern seal leakage is based upon the 
total gap (both initial gap and gap due to v.'aves) that is 
created at the stern seal. The amount of stern seal leakage 
is then dependent on the total gap and the plenum pressure 
that exists. It is recommended that the computer program 
be revised so that the amount of stern seal leakage is 
listed in two parts: leakage due to the initial gap and 

any additional leakage due to a gap being created by the 
waves . 



135 



^ . Steady-State Speed 

A minimum initial speed of 35 knots was used based 
on the results obtained in Pig. 2 and Appendix A. No 
computer runs were made to obtain the minimum critical 
speed between 30 and 35 knots under different vjave conditions, 
both ahead and astern. It is recommended that an investiga- 
tion into the minimum critical speed that can be used in 
the computer program be made to determine vjhat effect 
different wave conditions will have on the values of this 
minimum critical speed. 

5 . Validity of Computer Data 

Little or no performance data on the SES 100-B is 
available here at the Naval Postgraduate School. To check 
the validity of the computer data obtained in this study, 
a comparison of the performance data of the SES 100-B 
should be made under similar sea conditions as used in 
this study. A complete evaluation of the computer program 
for the selected variables that were studied could then be 
made . 

6 . Additional Data 

The range of vmve conditions was very limited in 
this study. Although basic trends could be developed, it 
is recommended that additional computer runs be made using 
simple waves from abeam and then using complex waves composed 
of 2 or more different waves. Realistic sea conditions 
should be modeled by the use of multiple simple v/aves as 
input data. 
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APPENDIX A 



CALCULATION OF CRITICAL SPEED 
From Nev7ton's Second Lav/ it is known that 

P = ma (3.1) 

The longitudinal force, F , is equal to 

F = Thrust - drag (3.2) 

Considering thrust as being constant, then we need only 
be concerned with those terms involving u to find the 
critical value of speed. Since u is composed of those 
terms that are u and u~ dependent, Eq. (3.1) and 
Eq. (3.2) can be written as 

mu = “(k^' u^ + k^ u”^'^^) + Thrust (3.3) 

Linearizing Eq. (3.3) by letting u = u^ + 6u gives 

-m 6 u = k^(u^^ + 2u^6u + 6u^) 

+ k^ u "1-56 _ Thrust (3.^) 

2 o 

Applying linear approximation, since we are interested 
in a steady-state critical speed for u^, and by specifying 
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that u ^’5^ r y , Eq. (3*^) can be v/ritten as 



-m6u = (2k^ u^) 6u + (1 + u^) (3-5) 



Letting X = — , then 



(1 + — ) 

u 

o 



-3/2 



(1 + x) 



172 



(3.6) 



Noting that 



(1 + x)^/^ [(1 + x)^(l + x)]^^^ 



(3.7) 



172 



(1 -I- x) (1 + x) 



and that the value of x, v/hlch Is equal to ^ , v;lll be 



o 



much less than 1, Eq. (3.7) can then be v:rltten as 



i tTo = =l-x (3.8) 

(1 + x)(l + x)'^ (1 + x)(l) 

Thus Eq. (3.5) can be written as 

• _q /p 

-m6u = (2k^ u^) 6u - 

= [(2k^u^) - 6u (3.9) 
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For 6u = 0, the critical speed, u , is from Eq. (3*9) 

o c 



2k^u = k^u 
1 oc 2 oc 



-5/2 



(3 



Solving for u gives 

O G 



k, 2/7 

%c = 



(3 



The values of k^ and k .2 are determined from the calm 
v/ater forces listed in Table III for 30 knots. 



k^ = FXPWAV = 5.898 X 10^ 

k^ = Total Forces - FXPWAV = 3.86? x 10^ 

Converting k^ and IC 2 into feet per second gives 

, _ 5.898 X 10 ^ 

2 (50..67)-3/2 

V = 3.867 X 10 ^ 

^ (50.67)^ 

where 50.67 feet/second = 30.0 knots. 

Solving Eq. (3.11) for gives 

u = 51.5 feet per second = 31.0 knots, 
oc 



. 10 ) 



. 11 ) 
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CONTROL STATEMENT 01400 
014 



CONTROL STATEMENT 01500. VALUES IN FLOATING POINT FORMAT. 

THIS CONTROL STATEMENT WAS NOT USED IN THIS STUDY SINCE STATEMENT IS CONCE 
WITH SHEAR AND MOMENT PROGRAM DATA, 
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